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We investigate the dielectrophoretic motion of charge-neutral, polarizable nanocolloids through molecular 
dynamics simulations. Comparison to analytical results derived for continuum systems shows that the discrete 
charge distributions on the nanocolloids have a significant impact on their coupling to the external field. Ag- 
gregation of nanocolloids leads to enhanced dielectrophoretic transport, provided that increase in the dipole 
moment upon aggregation can overcome the related increase in friction. The dimer orientation and the exact 
structure of the nanocolloid charge distribution are shown to be important in the enhanced transport. 
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I. INTRODUCTION 

Manipulation of microscopic particles with external electric 
fields, electrokinetics, has gained considerable attention due 
to applications in nanotechnology and biomedical research 
|H IjjIH 0]. One of the central methods is dielectrophoresis 
(DEP), the motion of polarizable particles due to a coupling 
with a non-uniform electric field J5[]. Recent studies have at- 
tested the versatility of DEP manipulation with carbon nan- 
otubes [6], DNA [7] and biological micro-organisms |8|], as 
well as the assembly of nanowires [9] and biosensor arrays 
JToh . The relevance of DEP methods to more complex inte- 
grated microscopic analysis and assembly devices 111 111 is evi- 
dent. 

It has been demonstrated that particles down to a few nm 
in size can be manipulated with DEP [.12, B, but it is still 
not known what is the minimum size of particles that can be 
efficiently transported and trapped |4|. Experiments are ham- 
pered by the fact that in situ tracking of nanoscale particles is 
very difficult. With decreasing particle size the effect of ther- 
mal noise becomes increasingly important. Furthermore, the 
motion of polarizable particles in non-uniform electric fields 
is inherently a non-equilibrium process, and particles in close 
contact may interact via non-additive many-body interactions. 
Given these challenges, it is fair to say that understanding the 
problem of nanoscale DEP calls for new theoretical insight. 

For an isolated particle much smaller than the characteristic 
length of the electric field non-uniformity, the DEP force is 
given by the well-known expression yfl 

Fdep = ip ■ V)E, (1) 

where p is the particle dipole moment induced by the field 
E. As the separation between particles decreases, their mu- 
tual interactions lead to modifications in their electric polar- 
izations and the DEP force affecting them J2, 0, EH] . It has 



'Electronic address: emppu.salonen@hut.fi 
^Electronic address: vattulai@csc.fi 
t Electronic address: mkarttu@uwo.ca 



been suggested that such changes result in enhanced electroki- 
netic trapping efficiencies at non-dilute conditions [ 16]. Con- 
trolled particle aggregation may in fact be an efficient way to 
overcome the large thermal forces that hinder DEP trapping 
and transport. Despite their importance, these issues, to our 
knowledge, have not yet been studied neither experimentally 
nor theoretically in detail. As for non-equilibrium simulations 
of DEP, there is only one recent study focusing on the trans- 
port of an individual nanocolloid 11711 . 

In this letter, we consider DEP through theory and model- 
ing and show how enhanced DEP transport can be obtained by 
controlled complexation. We use molecular dynamics simula- 
tions to study the aggregation and subsequent DEP transport 
of a dimer comprised of two spherical nanocolloids. We com- 
pare the transport properties of a dimer to those of a single 
nanocolloid, and analyze the competing roles of DEP force 
and friction. 



II. SIMULATION MODEL 

Our model lfl7tl is summarized as follows. We start from 
a charged spherical macroion with ATj small microions, each 
of charge q = +2e, electrostatically bound on the macroion. 
The macroion charge was set to Q = — A^g, so that the 
macroion-microion complex, termed hereafter as nanocol- 
loid, is charge-neutral. For excluded volume interactions we 
used a shifted repulsive Lennard-Jones potential, $wca = 
4e[(o-/{ry - r }) 12 - {oj{r VJ - r }) 6 + 1/4], where e and 
<7 are the characteristic energy and length of the interaction, 
respectively, and = \rj — is the center-to-center dis- 
tance between particles i and j. The potential was truncated at 
r« = r + 2 1 / 6 cr. The hard-core radius tq describing particle 
size was assigned a non-zero value in interactions involving 
macroions. Electrostatic interactions were calculated directly 
from Coulomb's law using a constant relative permittivity e r , 
given by (47reoe r ) _1 = 56 ea/e 2 , where erj is the vacuum per- 
mittivity. These types of generic models are physically trans- 
parent and have previously been employed in studies of elec- 
trokinetic phenomena in colloidal solutions [18, 19, 20, 21]. 
Our choice of simulation parameters can be viewed to model, 
e.g., DEP of reverse micelles [22], reverse phase emulsions 
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FIG. 1 : Dipole moments for two nanocolloids aligned parallel to the 
electric field, Eo — 0.96 e/ea, normalized with the dipole moment 
of a single nanocolloid. The symbols correspond to Ni = 10 (open 
squares), 100 (filled circles), and 300 (open diamonds). The solid 
line is the analytical continuum result 1 14], Inset: Total nanocolloid- 
nanocolloid attractive force for Ni = 10 (open squares), normalized 
with the value at D = 10 ao. The solid line shows the pure dipole- 
dipole contribution to the electrostatic force. 

JUt], or high-generation dendrimers [24]. 



HI. RESULTS AND DISCUSSION 

We first studied the influence of nanocolloid interactions 
on the dipole moments and the resulting electrostatic forces. 
The calculations were carried out in the absence of solvent. 
We considered the two specific cases where the line join- 
ing a pair of symmetric nanocolloids was either perpendic- 
ular or parallel to the direction of a constant electric field 
Eq = 0.96 e/ea. We used hard-core radii ro = 4.5, 16.9, 
and 29.6 a for the macroions, with a fixed microion surface 
density in all cases. The number of microions then increased 
with ro (Ni = 10, 100, and 300, respectively). The macroions 
were set a distance D apart from each other and subsequently 
held fixed at these positions. The microions were allowed to 
move freely, and their surface configurations were relaxed by 
slowly quenching the kinetic energy from the system with the 
Berendsen thermostat If25ll . 

Figure Q] shows the nanocolloid dipole moments for a pair 
aligned parallel to the electric field as a function of the 
nanocolloid separation D (in units of ao = ro + a/ 2). The 
solid line shows the analytical prediction in the continuum 
limit 11411 . For maximal polarization, we chose the value 6=1 
in eq. (7) of ref. 1 14]. For the largest macroion (Ni = 300) the 
enhancement in the dipole moments is in excellent agreement 
with the continuum case. Though, as expected, the agreement 
becomes worse with decreasing macroion size. This results 
from the small number of microions and pronounced discrete 
nature of the microion distribution. Yet it is remarkable how 
well the continuum theory grasps the essential features even 
just for Ni ~ 100. The results for the case of nanocolloids 



aligned perpendicular to the electric field also followed the 
same trend (data not shown): reduction of the dipole moment 
with decreasing interparticle distance, as follows from eq. (8) 
in ref. [14], was the most significant for Ni = 10. Again, 
for Ni = 300 our simulations were in an excellent agreement 
with the analytical result. 

The inset in fig.[T]shows the electrostatic force between the 
nanocolloids (Ni = 10) as a function of distance, together 
with just the bare dipole-dipole contribution. As the interpar- 
ticle distance decreases, the total electrostatic force reaches 
a value over an order of magnitude larger than the dipole- 
dipole contribution. By running a series of consecutive ki- 
netic energy quenching simulations at decreasing values of D, 
an electrostatically bound dimer was formed. Some of the mi- 
croions were trapped in the region between the macroions. For 
Eo > 0.96e/ecr, the resulting dimer dipole moments were a 
factor of 3 .4 - 4. 1 higher than for a single isolated nanocolloid. 
At Eo = 0.68 e/ea the ratio of dimer and single nanocolloid 
dipole moments was only 2.1. The above highlights the spe- 
cial nature of nanosized colloids: the smaller the charge, the 
stronger is the effect on creating asymmetries as the discrete 
nature of charge becomes important 112 111 . Here, as the mi- 
croions are strongly bound, that leads to non-uniform distri- 
bution that enhances the dipole for the pair. Note that a naive 
estimate, based on the notion that polarization is proportional 
to volume, would imply only an increase by a factor of two 
for the nanocolloid dimer. Here, at the higher values of Eq, 
the increase in the dipole moments is much larger. 

To fully describe the dynamics, we next simulated the DEP 
motion of a dimer with ro = 4.5 a and Ni — 10. We chose 
a spherically symmetric field, E(R) — Eq(Rq / R) 2 e^, where 
R is the radial distance from the electric field interaction cen- 
ter, Ro — 1500 a the characteristic length of the field, and 
e*R is the radial unit vector. To model DEP transport in a 
mesoscopic electric field geometry, with proper changes in 
the electric field strength, we employed an additional electric 
field system frame of reference and a periodic particle shifting 
scheme lfl7tl . The initial center-of-mass (CM) position of the 
dimer was set at R = Ro, with the line joining the macroions 
parallel to the electric field. A detailed description of our sim- 
ulation protocol is given elsewhere fl7il . 

The dimer was immersed in an explicit solvent of neutral 
particles with a density of p s — 0.3 a~ 3 . A cubic simulation 
cell of side length 45 a was used. Tests with larger systems 
showed that this cell size resulted in minimal finite-size ef- 
fects. The solvent density in our coarse-grained model is dif- 
ferent from the ones used in simulations of atomic or molecu- 
lar liquids, where typically p s w 0.8 a~ 3 . The solvent density 
in this letter was also used in the molecular dynamics study 
of electrophoresis by Tanaka and Grosberg 01811 . as well as in 
our previous paper on nanocolloid dielectrophoresis lfl7ll . The 
motivation for using p s = 0.3 a~ 3 , is that by reducing the sol- 
vent friction, related to p s , it is possible to study the effect of 
the spatial variance of the electric field to the DEP transport 
with simulations of reasonable length. Second, as we compare 
the DEP displacements of single nanocolloids and dimers [see 
below eq. ©], the solvent density has no major importance to 
the ratio of the displacements, and thus, to the main conclu- 



3 




0.5 1.0 1.5 2.0 

E (e/ea) 



FIG. 2: Mean radial DEP displacements at (a) T = 0.03, and (b) 
1.0 e/ks- The open (solid) squares show the data for dimers (single 
nanocolloids). The open circles at Eo = 0.96 e/ea show the data 
only for the simulations where the dimer assumed the higher dipole 
moment charge configuration (see text) for the whole length of the 
simulation. The solid lines show the predictions for the dimer case, 
see eq. 



sions of this letter. 

For each combination of Eq and temperature T we carried 
out 50-160 independent DEP simulations of length 5785 r, 
with the unit of time given by r = a^Jm/e. To give a con- 
crete idea of the timescales involved, with the central quanti- 
ties in our model having plausible nanocolloid values a = 0.2 
nm, m = 30 amu (atomic mass unit), e = fceTo, and To = 300 
K, our simulations would then correspond to a timescale of « 
4 ns. No dimer breakups were observed in the simulations. 

Figure |2] shows the dimer and single nanocolloid [26] CM 
mean radial displacements in the electric field system of ref- 
erence for the two extreme temperatures used, T = 0.03 
and LOe/fce- At Eq > 1.35e/e<7 the dimer displacements 
are clearly larger than those of the single nanocolloid. For 
Eo = 0.68 e/ea the DEP displacements are the same within 
the error bars. 

For Eq = 0.96 e/ea it was observed that the dimer polar- 
ization relaxed to two different states, corresponding to two 
different configurations of the microions trapped between the 
macroions. In some simulations the dimer assumed a state of 
very strong polarization with a dipole moment a factor of 5 
higher than that of a single nanocolloid. The fraction of sim- 
ulations with this strong dimer polarization increased with T 
from 15% at the lowest T to 33% at the highest T used. For 
the other, more common, polarization state the dipole moment 
ratio was only 2.1, i.e., a value similar as for Eq — 0.68e/etr. 



A more detailed analysis showed that in some cases the dipole 
made transitions between the two states. However, the dimers 
never retained the original polarized state that was formed in 
the absence of the solvent. Similar transitions were not ob- 
served for any other value of Eq. 

Equation ([1) shows that the DEP force affecting a dimer 
is larger than that on a single nanocolloid. But the frictional 
force due to the solvent on the dimer is also larger. Figure [2] 
shows that in our case the DEP force acting on the dimer out- 
weighs the increase in friction. Whether the DEP transport 
of aggregates is enhanced in comparison to single particles 
in general is essentially determined by the relative changes 
in these two factors. An estimate for the ratio of the friction 
factors of a dimer and a single spherical nanocolloid, £2/^1, 
can be made from hydrodynamics of smooth particles. With 
the long axis of the dimer aligned with the direction of mo- 



tion, £ 2 /£i ~ 1-29 B27I1 . For motion perpendicular to the long 
axis of the dimer, ^2/^1 ~ 1-43. This strongly suggests that 
nanocolloid aggregation enhances DEP transport, in line with 
our results in fig.|2j see also discussion below. 

To elaborate this issue, note that it has previously been ver- 
ified ifnll that DEP displacements with the range of external 
parameters Eq and T used here are linearly proportional to 
the mean DEP force. A simple estimate for the dimer DEP 
displacement, Ai?2, is then obtained from, 



AR 2 



F 2 (R ) /a 
Fi(Ro) \4 



Ai?i, 



(2) 



where Fi (Rq) is the initial DEP force on the system consisting 
of i nanocolloids at Rq, and Ai?j its displacement at the cor- 
responding Eq. The resulting values are shown as solid lines 
in fig. 12 The agreement with the simulation data at T = 0.03 
e/k-Q is good, though eq. (|2]) does not properly account for the 
spatial dependence of the DEP force. 

At increasing temperatures it should be expected that the 
rotational motion of the dimer hinders its DEP transport. This 
is due to two effects. First, for orientations deviating from full 
alignment with the direction of motion the effective friction 
factor of the dimer increases, see above. Second, the dipole 
moment of the rotating dimer changes as the microion distri- 
bution relaxes according to the field. Although changes in the 
microion surface distributions resulted, on the average, in ap- 
preciable DEP forces at all orientations, perpendicular dimer 
orientation with respect to the direction of E decreased the 
coupling strength up to by a factor of two. Considering the in- 
creased friction (by a factor of ~ 1 .43) together with the relax- 
ation of the microion distribution, decreasing the DEP force 
by 50%, an increase of -Fdep by a factor of ~3 is required 
for enhanced DEP transport of the nanocolloid dimer when 
rotations take place. 

Interestingly, the dimer serves as a model to gain insight 
into the role of T on rotational DEP motion of rod-like 
molecules. The external field counters the dimer rotation by 
exerting an electrostatic torque T = P2{R, 0) x E(R), where 
the dipole moment of the dimer, p 2 , depends on the radial po- 
sition R and the orientation angle with respect to the electric 
field. The characteristic time of dimer rotation may be shorter 
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FIG. 3: Electrostatic torque V on the dimer as a function of the 
dimer orientation, Eo = 1.91 e/ea. The solid line shows the results 
calculated for a relaxed charge distribution. 



torques are induced, aligning the macromolecules efficiently 
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FIG. 4: Distributions of the dimer orientations with respect to E in 
simulations with Eo = 1.35 e/ea and T = LOe/ks- Inset: values of 
the electrostatic torque V for different microion valences. 



than the charge relaxation time. This affects the actual values 
the dimer dipole moment assumes at different orientations. To 
illustrate this, we show in fig.[3]the values of torque affecting 
the dimer. Since the electric field is spatially variant, the data 
were collected only for R — 1464- 1539 a, where \E\ did not 
differ by more than 5 % from its value at i?o- We also calcu- 
lated the torque for fixed dimer orientations in the absence of 
solvent. The resulting charge distributions then corresponded 
to the case where the dimer rotation is very slow compared 
with the microion relaxation time. 

For low T the dimer rotations are slow and result in smaller 
absolute values of T. Yet, in comparison to k-gT, the values 
of the torque are extremely high and the dimer is efficiently 
aligned in the direction of the electric field. As T is increased, 
fluctuations in the microion distributions result in higher val- 
ues of T even for small 9. For lower values of Eq full rotations 
of the dimer were also observed. 

In addition to the case q = +2e and Ni = 10, we further ex- 
plored the effects of the microion valence by using two other 
microion configurations, namely q = +le, Ni = 20 and q = 
+3e, Ni = 7. The orientational distributions of the dimers are 
shown in fig. [4] Note that the distributions are biased toward 
cos 9=1 due to the initial orientation of the dimer in the sim- 
ulations. At E = 1 .35 e/ea and T = l.Oe/ks the DEP dis- 
placements for q = +le and +3e were 19% and 16% smaller, 
respectively, than for q = +2e. These differences were statis- 
tically significant with respect to the margins of error. The in- 
creased dimer rotation for the two extreme values of valence is 
also indicated by the smaller values of the electrostatic torque 
affecting the dimer (cf. inset of fig. [4J. These results demon- 
strate that the microion valence could in fact be a control pa- 
rameter for nanocolloid DEP motion. 

So far rotational effects under DEP have been considered 
only for much larger objects with large aspect ratios, such 
as carbon nanotubes [6] and DNA 1 28] . In these cases large 



along the electric field. However, in our case even small per- 
turbations of the microscopic charge distribution due to ther- 
mal fluctuations can result in considerable changes in T(R, 9), 
see fig. [3] The rotations could lead to non-negligible hydro- 
dynamic coupling between DEP-manipulated dimers or larger 
aggregates. Though, as previously mentioned, our modeling 
can be viewed to correspond to a timescale of a few ns. Since 
we model a dilute system, possible hydrodynamic effects be- 
tween rotating dimers are thus not relevant here. 



IV. CONCLUSIONS 

We have shown how the proximity of nanocolloids with dis- 
crete charge distributions increases their DEP coupling. By 
forming nanocolloid dimers, enhanced transport can be ob- 
tained provided that the increase in the dipole moment over- 
comes the increase in the friction. We have found the rota- 
tional motion of dimers to hinder the transport due to smaller 
DEP coupling and larger friction. The microion valence was 
shown to affect dimer orientation and transport. Our results 
indicate that controlled aggregation can be an efficient way of 
overcoming thermal forces and obtaining high-precision DEP 
particle manipulation at the nanoscale. 
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